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ABSTRACT 



A situation is considered in which a grenade launcher fires at a target, 
observes the projectile impact, adjusts the point of aim, fires a second 
round, adjusts, fires, and so on. A model is developed, based on the 
assumptions of normal ballistics, perfect observation of impacts, adjust- 
ments without error, and a unimodal target destruction function. The 
problem is to determine optimal adjustments, in order to maximize the 
probability of target destruction within a given number of rounds. It 
is shown that seemingly different adjustment procedures are equivalent, 
if viewed in appropriate coordinate systems. Previous results concerning 
sequential adjustments which are constrained to be linear functions of 
observed impact points are extended to the class of translation invariant 
procedures. Properties of the optimal sequential adjustment procedure, 
including some related to stochastic approximation, are reviewed. The 
effects of errors in judging impact positions are discussed. 
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1. Introduction, 



Consider the problem of adjusting grenade launcher fire on a target, 
using a "shoot-adjust-shoot” strategy. After each round is fired, the 
grenadier observes the projectile impact, and estimates its distance from 
the target. He then attempts to adjust his point of aim, so as to place 
the mean of the ballistic distribution of impact points directly over the 
target. The problem we address is that of determining an optimal adjust- 
ment scheme. 

In what follows, we assume the ballistic distribution of impacts 
of rounds fired at a given aimpoint is bivariate normal. We further 
assume that the target is a point target, and that the target destruction 
function (or lethality function) d is the same for all rounds fired. 

The value of d(x,y) is defined to be the probability that the target 
is destroyed by a round impacting at the point (x,y) relative to the 
target. The difficulty of finding optimal adjustment procedures hinges 
initially upon how errors are made by the grenadier in estimating the 
position of each impact point, relative to the target. For the most part, 
we confine ourselves to the simplest case in which the grenadier is 
"perfect” — he delivers adjustments which are exactly the desired propor- 
tions of the observed miss distances. In the last section, some remarks 
are made concerning the probable effects of errors in this process. 

The problem of how to optimally adjust fire on a stationary point 
target has been of considerable interest. Grubbs [8] gives an adjustment 
scheme which is "best” (i.e. minimum variance) and unbiased among linear 
functions of the impact points. (In the literature of statistical 
estimation, such a scheme is called "BLUE”, for "best JLinear unbiased 
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estimator. M ) Nadler and Eilbott [10] derive the BLUE adjustment scheme 
and show it maximizes, among linear functions of the impact points, the 
probability of destroying the target with each round, for a fairly general 
class of destruction functions. This procedure is a special case of 
stochastic approximation schemes, about which a great deal, mostly of an 
asymptotic nature, has been written [5], [11], [12]. Other stochastic 
approximation approaches to the adjustment of fire have included cases 
in which the observer cannot precisely locate the impact points (see, 
for example, Barr and Piper [3]). 

The main purpose of this paper is to demonstrate that, with appro- 
priate coordinate systems, the procedures of Grubbs and of Nadler and 
Eilbott may be obtained from well known statistical theory concerning 
BLUE estimators. We also show (in Section 2) that these procedures are 
indeed stochastic approximation procedures, and (in Section 3) that they 
possess an optimality property stronger than BLUE. The present approach 
seems notationally and conceptually simpler than those used previously. 

It is notationally convenient to restrict our consideration 
to a single dimension, say range, in what follows. Under the assumptions 
we make below (which are somewhat weaker than those of Nadler and Eilbott) , 
the restriction to one dimension causes no loss of generality: identical 

optimal one-dimensional procedures carried out Independently in each 
dimension gives an optimal two-dimensional procedure. (See Section 5 below.) 

2. The Model . 

Assumptions which permit analysis of a two-dimensional problem by 
separate analyses with each dimension are essentially those required for 
Anderson’s [2] inequality. The following are sufficient: 
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i) statistical independence of the range and deviation components 
of the impact points, each component having a symmetric uni- 
modal continuous distribution, 

ii) no grenadier errors are made, and 

iii) the destruction function d is symmetric about the impact point 
(x,y) and {(x,y): d(x,y) > a} is convex for all a ^ 0. 

In particular, a circular normal ballistic distribution with a 
destruction function of the "circular normal" or "circular cookie-cutter" 
types will satisfy condition (iii) . It is well known that in actual 
gunnery problems, the ballistic distribution is such that the variance 
in range is not the same as that in deviation. In most gunnery problems, 
these variances depend on many factors, but principally upon the grenade 
launcher-target range. The assumption that the impact points fall in a 
circular normal pattern about the aimpoint is equivalent to agreeing to 
measure range and deviation in the corresponding standard deviation units, 
and in addition, assuming independence between the range and deviation 
components of the impact points. 

Suppose, then, that X^,X^,X^,... are independent normal random 
variables with unknown mean y and known variance a 2 (we may take 

9 

o^=l without loss of generality ) . These random variables represent 
the sequence of impact point ranges that would be obtained with a sequence 
of rounds fired without adjustment (i.e. with a fixed aim point). The 
problem is to determine a sequence of adjustments, m * m 9 w ^ ere 

is the adjustment made between the and n + 1— rounds, so that 



As mentioned above, assuming a 2 38 1 is equivalent to measuring distances 
in a-units. In theory, a 2 varies with the range being fired, and hence 
varies as adjustments in range are made. However, in practice, this 
variation in a 2 is sufficiently small within a given firing mission that 
a 2 may be assumed to be constant. 
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the probability of target destruction with each round is maximized. (As 

we remark below, this will maximize the probability of target destruction 

within any fixed number of rounds.) 

t h 

In general, the n — adjustment is a function of the first n 

observed impact ranges, X^,X^, . . . ,X^. Each adjustment £ has the 

effect of translating the coordinate system units. For example, if 

an adjustment £ = -50 (corresponding to the intention of dropping 50 

meters in range) is made, the net effect is as if the coordinate system 

(i.e. target) were transformed 50 units away from the grenade launcher. 

We must therefore keep track of the adjustments made up to the n-^- 

th 



round, so as to be able to express the 



impact point in coordinates 



which are meaningfully related to the previous impact points. Accordingly, 
let X^ denote the j— impact range, expressed in the r^* coordinate 
system. Upon firing the first round, the grenadier observes X*, and 
makes adjustment The first round can thus be given as X^ = x| + 

Similarly, the second round X^, fired after the adjustment 5^, can be 
given in the initial coordinate system as X* = X^ - and so on for 

subsequent rounds, adjustments and coordinate systems. 

Remark . It should be noted that the procedure we will determine 
does not require the grenadier to actually estimate the miss distances X^. 
Rather, it will turn out that he should attempt to adjust by a fraction 
of the observed miss distance. Thus, the grenadier needs only to be able 
to duplicate some fraction of an observed distance in his aim adjustment. 
The present discussion involving estimates X^ of actual distance values 
are a theoretical convenience only. 
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If all observed impact ranges are transformed to the initial 

coordinate system, and if each adjustment is some linear function of the 

data available for that adjustment (i.e. . ,X^) ) , then the 

5*s are also linear functions of the "unadjusted 11 impact ranges 

xj ,xj; , . . . jX 1 . The problem then becomes that of finding the optimal 
1 L n 

linear function of a random sample of identical, independent normal 
(y,l) random variables. Naturally, one expects the sample mean 

1 n ] 

X 1 = — I x ] , to provide an optimal adjustment -E 1 , from the initial 
n n 1 n 

aimpoint y , where we write "5^," rather than to indicate we 

mean adjustment from the initial aimpoint. In what follows, we call this 
the "X-bar procedure," 

It is easy to establish that, in fact, = -X^ is the optimal 

linear adjustment. It is well known that X^ is a BLUE estimator for 

* 

y under the present conditions. With the unimodal, symmetric target 

destruction function considered here it is clear that one may restrict 

his attention to linear adjustments which are unbiased estimators of -p, 

since otherwise the probability of target destruction on the next round 

could be made as high or higher by subtracting the bias. But, among the 

s t 

normally distributed aimpoints for the n + 1 — round, which arise after 
applying unbiased adjustments the probability of target destruction 

is maximized by minimizing the variance of the aimpoints about the target. 
That is precisely what the BLUE estimator of -y achieves, and it is the 
optimal scheme developed by the aforementioned authors. 

For example, Larson [9], gives the following [Theorem 7.3.1] s X is a 
random variable with mean y and variance a^. If a 

random sample of X, X is the best linear unbiased estimator of y. 
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In order to see that the X-bar procedure {E 1 } is a stochastic 

n 

approximation procedure, we shall express it in terms of the sequence of 
current coordinate systems, Note that 



_y 1 = _vl 

X 1 V 



After making this adjustment and firing the second round, resulting in 
X|, we have, in the second coordinate system, two observations: X^ and 

X| . Now 

£■2 = rl _ r 1 

t, 2 ^2 H 

- - \ (Xj+xj) - 

= ” 2 



In general, 

= - ^ (Xl+(n-l)^ ,) - , 

n n n-1 n-1 

- - i x". 

n n 

Thus, the optimal M X-bar n procedure is identical to the M — ■ - 

stochastic approximation" [11] procedure which, after the nr^- round, 

adjusts by the fraction of the observed miss distance of the n^ 

round. The good asymptotic properties of such stochastic approximation 

procedures are not surprising when viewed in this light. Similarly, 

the distribution of the n + 1— impact point derived in [10] 

n “t i 



using multivariate analysis, is now very easily seen to be normal with 
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mean zero and variance (1 H )o 



U 9 n + 1 
_) 0 z . 

n n 



. That the mean is zero follows 



from the fact that X is unbiased for y = E(X* ). The variance of 

n rrt-l 

the aimpoint for the n + 1— round is V(X^) = anc ^ since X * | ^ is 



independent of X* V(x£jj) = VCX^) + V(X*) - 1 + K 

Remark . It is easily shown, by induction for example, that the 
n— adjustment in terms of the j— coordinate system (data and shift) 

1 i n 

Is that which places at zero (i.e. on the target); £ * ~ *T I X . . 

n n u - L 

The earlier observation that £ n = - — X n also follows from this, since 
„ n n n 

n-1 

I X* = 0. 

i=l 



3 . Strong Optimality of the X-bar Procedure . 

We have re-established, by elementary methods, that the X-bar 
procedure has, under our assumptions of independence, normality and 
regularity of the target destruction function, the properties 

(a) among linear adjustments, it maximizes the probability of 
destroying the target with each round fired, and 

(b) it provides unbiased aimpoint s (after the initial round is fired) 
with minimum variance. It easily follows that also it 

(c) maximizes the probability of destroying the target within the 
first k rounds; k = 1,2,..., and 

(d) minimizes the expected number of rounds required to destroy the 
target. The latter may be established using the identity [4] 

QO 

E(N) = l P[N * n], 

j-1 

where N is the (random) number of rounds required to destroy 
the target.Evans [6] gives an elegant proof of the intuitively 
obvious fact that maximizing the probabilities of target destruc- 
tion with each round minimizes the expected number of rounds 
required to destroy the target. 

Since 'x is a robust estimator of y, with central limit theorem 
properties when the X^s are not necessarily independent and identical, it 
follows that the X-bar adjustment also, in loose terms 
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(e) is robust (i.e. tends to operate efficiently even when the 
underlying assumptions are not met) , and 

(f) tends to give normally distributed aim points. 

We wish to show next that the X-bar procedure is optimal in a 
class of procedures much wider than the class of linear procedures 
considered heretofore. The approach is to use a result of Wolfowitz [3], 
concerning minimax estimators, in a decision theoretic setting. We use 
the notation and terminology of Ferguson [7]. In the preceding section, 
it was established that the problem of finding an optimal aimpoint after 
n rounds have been fired is equivalent to finding an optimal estimator 
for y, based on a random sample from a N(y,l) popula- 

tion; once that estimator y is found, the optimal adjustment for the 
next round is -y , from the initial aimpoint. Clearly, this problem is 

* A 

invariant under the group of translations (i.e. roughly speaking, if y 

is optimal for y, y + c is optimal for y + c) . It is, therefore, 

extremely reasonable to insist that the solutions to the problem should 

be invariant (i.e. roughly speaking, if y(xJ,...,X 1 ) is optimal for y, 

1 n 

then y (X}+c,Xj:+c, . • • ,X*+c) should be optimal for y + c) . We shall, 

1 l n 

therefore, restrict our attention to estimators which are invariant under 
translations; this means essentially that we are indifferent as to where 
the origin of our coordinate system is placed, relative to the target, 
before firing begins. (It is intuitively unacceptable for an arbitrary 
assignment of this sort to affect the procedure to be followed in adjust- 
ing fire.) 

'k 

For a discussion of invariant problems and solutions, see Ferguson [7], 
Definitions 4.1.2 and 4.2.1. It should be noted that, while functional 
equations of the form f (x+c) = f(x) + c determine linear functions in 
the single variable case, this is certainly not the case for functions 
of several variables (see Aczdl [1]). 
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Theorem . In the class I of translation in variant procedures, the 
X-bar procedure is uniformly best, in the sense that, for any t > 0, 

max P[|y(x} X*) - y| < t] 

Cei 1 

- P[|? - u| < t]. 

Our proof makes use of results of Wolfowitz, which are stated in 
convenient form by Ferguson (Theorem A .6.1) as follows: 

Lemma. If x} , . . . ,X* is a sample from a normal distribution with mean y 
1 n 

and known variance [1, say], then X^ is a best invariant estimate of y 
and a minimax estimate of y, provided that the loss function L(y,y) 
is a nondecreasing function of |y - y| and that E[L(0,X^)|y = 0] 
exists and is finite. 

To say that X is a best invariant estimate of y means that 
n 

it is as good as any other invariant estimate: for every y and any 

__ 

M € I , 

E L(y ,5c 1 ) s: E L(y,y(x 1 ,...,xh), 
y n y in 

where "E^ 11 denotes expectation with respect to the joint normal distri- 
bution in which the x} have mean y. To say that X 1 is minimax 

i n 

estimate of y means that, roughly speaking, it minimizes the maximum 
loss we could sustain, regardless of ’’natures choice” of y. 

Now let 

/ o if I y - y I < t 
L(y,y) = < 

' 1 otherwise. 
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Clearly L satisfies the conditions of the lemma, since for any fixed 
t > 0, 

i) L is nondecreasing in |p - p|> and 

ii) E L (0 ,X 1 ) = P [ | X 1 1 t|w = 0] £ 1. 

U n n 

By the lemma, It follows that, for any y, 

E y L(p,p) - P[|p - p| k t] 

is maximized over I by taking y = X^, which completes the proof. 

The theorem may be interpreted, along the lines of the argument 

in Section 2, as guaranteeing that the X-bar procedure 

(g) has properties (a) through (d) in the class of translation 
Invariant procedures . 

As mentioned above, it Is difficult to imagine that one would ever want to 
consider adjustment procedures outside the class I. It should be noted 
that the class I contains the class of linear procedures as a proper 
subset; indeed, intuitively speaking, the linear procedures form a "small" 
subset of I . Thus our theorem may be described in broad terms as dropping 
the "linear" restriction imposed by Grubbs and Nadler and Eilbott. 

Our application of Wolfowitz*s results can be extended to a larger 
class of target destruction functions, by taking the loss function to be 

L(y,y) = -P [Procedure using 5 1 = y(xJ,...,X ) 
y n 1 n 

destroys the target on the n+lst round] 

r 

= -d (p-p (X^ , . . . »X^)) 
l n 

The X-bar procedure will be best (in I), provided that the destruction 
function is a nonincreasing function of |y - y | ; existence and finiteness 
of E-d(-p(x} , . . . .X 1 )) is automatic. (Relationships between d(x,y) and 
d(x) are discussed below, in Section 5.) 



11 



4 . Summary and Conclusions , 

The foregoing analyses indicate that an optimal (in several senses) 
adjustment scheme for the grenadier to pursue consists of shifting the 
aimpoint by an appropriate fraction of the miss location of the preceding 
grenade burst. It is significant that this procedure does not require 
numerical evaluation of miss distances, nor does it require computations 
with, or manipulations of the sequence of impact points observed as firing 
proceeds. The procedure we have studied is well known and often used in 
practice — our contribution has been to determine its optimality under 
simplifying assumptions. 

The most stringent of the assumptions made in our analysis is that 
the grenadier can produce a shift in aim which is precisely some fraction 
of an observed miss distance. Some insight into the effects of relaxing 
these assumptions can be gained by considering the analytical expression 
in our model for the variance of aimpoint after n rounds. It can be 
argued that the grenadiers inability to exactly produce the desired shifts 
in aimpoint will have the effect of increasing variance in aimpoint. 

If the grenadier makes normally distributed "errors" with mean 0 
and variance t 2 , independent of the actual miss distance, the overall 
situation is as if the ballistic distribution has variance a 2 + t 2 = 1 + i 2 , 
with an errorless grenadier. That is, under the independence assumption, 
the mathematical model reduces to the previous case with a larger variance. 
The variance of the aimpoint for the n + 1— round is now (l+t)/n, 
as compared to 1/n with the errorless grenadier. 

A more realistic model for grenadier error is that the variance 
decreases with miss distance. Such would be the case, for example, if 
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the conditional distribution of grenadier error, gives a miss distance 
x, is uniformly distributed over the interval x + px, where p repre- 
sents a proportion parameter measuring the grenadier’s skill. This 
situation places a premium on obtaining rounds near the target, as soon 
as possible. It is anticipated that the X-bar procedure will be shown 
to be optimal with a variety of grenadier error models of this general 
type; we are continuing with work along these lines. Simulations with 
models of this type appear to support a claim that the variance in aim- 
point does not go down like 1/n, as it did in the earlier models we 

considered. Rather, the variance appears to go down somewhat slower with 

- 3/4 

the latter model; perhaps as n or so. 

These variations in the variance in aimpoint under changes in the 
grenadier error models and parameters should provide a method of assessing 
the value of multiple shot grenade launchers. This might be accomplished 
by assessing the benefits of the multiple round capability (vice a single 
shot capability) in terms of the grenadier error characteristics, which 
in turn affects the variance in aimpoint. The aimpoint distribution should 
still be approximately normally distributed, by the central limit theorem. 
Hence, with a knowledge of the aimpoint variance and ballistic variance, 
together with a target destruction function, one can compute the proba- 
bility of destruction with each round, as well as the expected number of 
rounds required to destroy the target. 
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5 . Appendix . 

We have asserted that an optimal two-dimensional procedure is 
obtained using optimal one-dimensional procedures in each dimension. 

The purpose of this appendix is to verify that this is the case, under 
the following assumptions: 

(i) The range X and deflection Y components of the impact point 
of a round fired at aimpoint (0,0) are independent, symmetric unimodal 
continuous random variables, with densities f (x) and g(y), respectively. 

(ii) The target destruction function d(x,y) symmetric in each argument, 
where d(x,y) denotes the probability that a round impacting at (x,y) 
destroys a target at (0,0). Equivalently, under these symmetry condi- 
tions, d(x,y) may be viewed as the probability that a round impacting 
at (0,0) destroys a target at (x,y). 

(iii) {(x,y): d(x,y) > a} is convex, for all a. 

The approach we take is to present two "auxiliary" destruction 
functions, related to d(x,y), which are one dimensional and satisfy 
the conditions of Anderson’s inequality. It is 9hown that, at each stage 
of the sequential adjustment procedure, optimal adjustment for the 
"marginal" problems will give the optimal "joint" adjustment. Suppose 
the target is at (0,0). The probability of destroying the target, given 
one aims at (y^,^) i s 



p (Pf ,^ 2 ^ 55 P [Target destroyed | aim at (v^,^)] = 



d (x,y) f (x+y 1 )g (y+u 2 )dxdy 



d (x-v^ , y-u 2 ) f (x) g (y) dxdy , 



_ 00 __oo 
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where f (x+y^)g(y+y^) is the bivariate distribution of impact points 
with aimpoint (y^,y 2 ). problem is to show that this is maximized 

with y^ = y 2 = 0; Nadler and Eilbott do so [Lemma 2] using Anderson’s 
inequality, which requires that d(x,y) be symmetric and {(x,y): d(x,y) > a} 
be convex for all a. Now 



00 

p(y 1 ,y 2 ) = 

oo 




d(x-p ,y-p )f (x)dx 



g(y)dy 



is maximized if it is possible to maximize 



d 2 (y-p 2 ) 



d(x-p 1# y-p 2 )f (x)dx 



over y^ for each value of y - y 2> independent of the value of y - y^, 
and then to maximize 



d 2 (y-y 2 )g(y)dy 



over y ^ • Let d^(x) = d(x,y). As a function of x alone, note that 
d^(-x) = d(-x,y) = d(x,y) = d^(x) ; and {x:d^(x) > a} is convex, for 
otherwise { (x,y) :d (x,y) > a} would not be convex. Thus Anderson’s 
inequality implies that, for the one dimensional problem with ballistic 
distribution f(x) and any of the destruction functions d^(x) , the 
maximum of d^ is attained with y^ = 0, independent of y^* Now apply 
Anderson’s inequality again for the other dimension, with ballistic dis- 
tribution g(y) and destruction function d^(y) . [Note that 



d 



2 




d(x,-y)f (x)dx 



d(x,y)f (x)dx = d 2 (y) , 



00 
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and {y:d 2 (y) > a} = {y 



d(x,y)f(x)dx > a} is convex, again because 
{ (x,y) :d(x,y) > a} is convex for all a.] It then follows that 



J d(y-y 2 )g(y)dy is maximized with y^ « 0, 

Thus, given any two dimensional problem satisfying the conditions 
posed above, an optimal solution is obtained following optimal one 
dimensional aim strategies at each stage. Note that radial symmetry of 
the bivariate destruction function d( # ,*) is not required. 
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